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ABSTRACT 

Binary radio pulsars are generally believed to have been spun up to millisecond periods 
i-C ■ (i.e. recycling) via mass accretion from their donor stars, and they are the descendants 

of neutron star low-mass X-ray binaries. However, some studies indicate that the 
q , formation of pulsars from the accretion-induced collapse (AIC) of accreting white 

!h ■ dwarfs (WDs) cannot be excluded. In this work, we use a population synthesis code to 

r/j ' examine if the AIC channel can produce eccentric binary millisecond pulsars (BMSPs) 

in the Galaxy. Our simulated results indicate that, only when the natal MSPs receive a 
relatively strong kick (> 100km s _1 ), can the AIC channel produce ~ 10 — 180 eccentric 
\ (e > 0.1) BMSPs in the Galaxy, most of which are accompanied by a Helium star. 

Such a kick seems to be highly unlikely in the conventional AIC process, hence the 
' probability of forming eccentric BMSPs via the AIC channel can be ruled out. Even if a 

£T) , high kick is allowed, the AIC channel cannot produce eccentric BMSPs with an orbital 

period of > 20 days. Therefore, we propose that the peculiar BMSP PSR J1903+0327 
' cannot be formed by the AIC channel. However, the AIC evolutionary channel may 

qq , produce some fraction of isolated millisecond pulsars, and even sub- millisecond pulsars 

■ if they really exist. 

Key words: binaries: close - pulsars: general - stars: white dwarfs - Galaxy: stellar 
content 
> ■ 

C3 ; 1 INTRODUCTION B When the mass transfer ceases, a binary millisecond pul- 

_ , , nrnn ^ , . , , sar (BMSP) is formed. 

Millisecond pulsars (MbPs) are characterized by snort pulse TT ,, . . „ ,, ^-n 

V ' ■ _ lq However, the origin ot the MbPs still presents some 

periods (P < 20 ms). low spin-down rates (P ~ 10 — . , , ^. , ^ TO rovn 

r _ 21 -i^ n ^ ■ ■ i P/(2P) controversial puzzles, lirstly, can the known Nb LMXBs 

a SS m n° °i arac erls ^ a S es \ T A ) evolve into the observed M SPs both in the Galactic disk 

10 — 10 yr , and weak surface magnetic fields (B ~ i T , „ . „ AT — — | | lr , 00 Ko .. , , . , 

„ q , r— r ; 1 I l ° „ \ , Kulkarm & Naravan 1988 ; btatistical analysises show 

10 8 - 10 9 G (Manchester 2004). About 80% of the * , , , , , , \ t hiyd ■ i o a t •+ a 

' M ■ 1 , — *• m r , , that the birthrate of LMXBs is ~ 1 — 2 order of m agnitudes 

Lf 6 m u bmar f s ' bu * on % n, 1% ' o , the total pu l- lQwer than that Qf Mspg J CotefcPvlvse J|i93a EoTinteil 

f rs dBhattacharya fc van den Heuvel || l99l| ;| Lorimeri | 2008| ) . ^ SecQndl - t ig hard tQ understand how isolated MS p g 

btandard models proposed that MbPs are formed m neu- . „ , ,. ,. , „ , . , , 

, ,r^,vr, \ n n m *h e Galactic disk were formed via the standard recycling 

tron star (Nb) low-mass X-ray binaries (LMXBs), m which . , ,. r ,, , n 

, ' J , ' scenario. Although evaporation ot the donor stars by the 

the Nb has accreted the mass and angular momentum from , . , ,. ,. „ , , m-o-n ■ 1 ^ a r n 

° nigh energy radiation ot the MbPs might account tor the 
its companion, and been spun up to a short spin period . . „ . , , , „ Dn , , . , r 

r- iH — -i I f- - 1 - ii — -'tt origin ot the isolated MbPs, the observed time scales for 

dAlpar et al. 1 1 19821 : iTauris fc van den Heuvell 120061 ). Mean- ,. . ' , ,, 

1 , ■ — — : ' evaporating a companion star seem too long, thirdly, m 

while, mass accretion onto the Nb induces its initial mag- , , , , , , . . , , „ . ,. . 

, „ ,„ ,„ „. „ o q „ standard model the strong tidal effects operating in the bi- 

netic field (B ~ 10 12 - 10 13 G) to decline to B ~ 10 —10 G 5 6 
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1 As an alternative mechanism of accretion-induced field decay, 
the magnetic fiel d of the NS may be screened or buried by the ac- 



crete d material l lTaam 
l200lf) . 



van den He uvel 1986; Cummin g et al.l 
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nary during the mass transfer serves to circularize the orbit. 
However, the discovery of the eccentric (e = 0.44) BMSP 
PSR J1903+0327 in the Galactic plane h as challenged the 
standard scenario jChampion et al. II2008T ). 

As an alternative evolutionary channel, BMSP may 
originate from the accretion -induced col l apse (AIC) of an 
accret ing white dwarf (WD) (|MichelHl987n .ll vanova fe Taaml 
|2004l ) suggested that AIC may occur by thermal time- 
scale mass transfer in such binaries with orbital periods 
of a few days. Once the accreting ONe WD grows to the 
Chandrasekhar limit, the electron capture process may in- 
duce gravitational collapse rather than type la supernova. 
As a result of angular momentum conservation and mag- 
netic flux conservation, an MSP with rapid spin and low 
magnetic field may be produced. To account for its observed 
characteristics, the bursting pulsar GRO J1744-28 was sug- 
gested to be originated from the AIC of a massive ONe WD 
(van Paradiis et al.lfl997l ; IXu fc Li 1120091 '). Recently, the es- 
timated birthrates of MSPs by population synthesis calcula- 
tions sh ow that the often-n eglected AIC channel cannot be 
ignored (|Hurlev et alj|20ich . 

In the AIC channel the puzzles in the standard re- 
cycling model mentioned above might disappear. Firstly, 
the AIC of accreting WD s has been raised to int erpret 
the birthrate discrepancy (|Bailvn fc Grindlav I [l990.) , and 
the AIC process may be associated with the r-process nu- 
cleosynthesis of heavy (ba ryon number A > 130) nuclei 
<|Qian fc Wasserburg (20031 ). Secondly, a kick velocity caused 
by asymmetric collapse may be produced during AIC of 
WDs. An appropriate kick can disrupt the binary system, 
and results in the birth of isolated MSPs. Otherwise, the 
bin ary survives and an eccentric BMSP is formed. 

iFerrario fc Wickramasinghd ([2007) argued that the 
AIC channel can form BMSPs of all of the observed types. 
IChampion et al.~l (|2008h proposed that the AIC of a massive 
and rapidly rotating WD could produce the observed orbital 
parameters of PSR J1903+0327. In this work, employing the 
binary population synthesis approach we attempt to inves- 
tigate if the AIC evolutionary channel can produce a popu- 
lation of eccentric BMSPs, as well as isolated MSPs. In Sec- 
tion 2, we describe the population synthesis approach and 
the evolution model of MSPs. The simulated results by pop- 
ulation synthesis are given in Section 3. Finally, we present 
a brief discussion and summary in Section 4. 



2 INPUT PHYSICS 
2.1 Population synthesis 

Using an evolutionary population synthesis based on the 
rapid binar y star evolution (BSE) code (e.g. iHurlev et all 
2000, 2002), we attempt to study the statistical properties 
(such as the birth rate, total number, the distributions of or- 
bital period and eccentricity) of MSPs formed via the AIC 
process of a massive ONe WD in the Galaxy. In calculation, 
we consider the evolution of single stars with binary-star in- 
teractions, which includes the mass transfer and accretion 
via stellar winds and Roche lobe overflow, common enve- 
lope (CE) evolution, supernovae and AIC kick, tidal friction, 
and orbital angular momentum loss containing gravitational 
wave radiation and magnetic braking. 



All stars are assumed to born in binary systems, and 
start with zero eccentricities and a solar metallicity (Z — 
0.02). We adopt the following input parameters for the sim- 
ulation. (1) A constant star formation rate S = 7.6085 yr _1 , 
which corresponds to that one binary with Mi ^ 0.8 M© is 
born in the Galaxy per year; (2) The primary mass distribu- 
tion $(lnA/i) = Mi£(Mi), in which the initial mass function 
£(Mi) is given bv lKroupa et all (|l993h : (3) The secondary 
mass distribution <3?(lnAf2) = M2/M1 = q, which corre- 
sponds to a uniform distribution of the mass ratio q = — 1; 
(4) A uniform distribution of lna for the binary separation 
a, namely 3>(lna) = fc = constant (no rmalization results in 
k = 0.12328) (see IHurlev et alj|2002h . The input parame- 
ter space for the primary mass Mi, the secondary mass M2, 
and the separation a are set to be 0.8 — 80Mq, 0.1 — 80M©, 
and 3 — 10000-Rq , respectively. For each initial parameters 
X (Mi,M2, and a), we set n x grid points in logarithmically 
space, so and have 

lnXmax — lnXmin 



<51nx 



(1) 



During the evolution, if a binary appears as BMSP, it makes 
a contribution to the birthrate as 



Sr = S$(lnMi)$(lnM 2 )$(lna)51nMi51nM2(51na. 



(2) 



If this system lives for a time At as a member of BMSPs, it 
makes a contribution to the number 



Sn = 5r x St. 



(3) 



During the CE evolution, the parameter qce = 
-EWd/(-Eorb.f — £?orb,i) describes the efficiency that the 
orbital energy is transferred to expel the envelope 
l|Hurlev et alJl2002h . here -Ebind is the total binding energy 
of the envelope, E OI h,i, and i?orb,i are the final and the 
initial orbital energy of the cores, respectively. There ex- 
ist so me diversity in calc ulating the initial orbital energy, 
hence IHurlev et all (|2010h suggest that in this formulation 
Qce = 3 s h ould be consistent with oce = 1 adopted by 
iPfahl et al.l (|2003h . Both the estimated velocity of individ- 
ual pulsars and the pulsar velocity distribution indicate that 
the nascent NSs should have received a kick velocity, which 
may ro ot in the asymmetric collapse of massive proge nitors 
of NSs (|Burrows fc Haves lll996l : lLai fc GoldreichhoOOl ). The 
kick plays an important role in determining the fate of bina- 
ries, whether surviving or disruptive. Noticeably, it is widely 
accepted that the NS formed by the AIC may receive a lower 
kick velocity (e.g. <taic = 20km s _1 in the standard model of 
Hurley et al. 2010). However, some works presented the op- 
posit e point of view l|Harrison fc Tademarulll975l ; IE ai et al.l 
120011 ) (see section 4). In our standard model, a relatively 
stronger kick dispersion for NSs by core collapse sup ernovae 
is taken to be a C c = 265km s" 1 i|Hobbs et al. |2005| )FI. while 
(Taic = 150km s _1 (a ultra-high kick diepersion) for NSs via 
the AIC. To study the influence of AIC kick on the formation 
of MSPs, other kick dispersion such as cr aic = 100, 50km s _1 
are also included. Unless we particularly mention , the BSE 
input parameters in table 3 of lHurlev etafl ([2002) are used. 



2 To match the observed intrinsic ratio between disrupted recy- 
cled pulsars and double NSs, the simulated results by population 
synthesis show that the natal kick of NS formed in c lose interact- 
ing binaries is ~ 170km s" _1 jBelczvnski et al.ll201ol) . 
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2.2 Evolution of MSPs 

When the mass of the accreting WD reaches the Chan- 
drasekhar limit, the process of electron capture may dom- 
inate, and induce the WD t o collapse to be an NS or ex- 
plode as a type la supernova l|Canal et al. Ill980l ; llsern et al.l 



Table 1. Model parameters for binary population synthesis. 



1983 ), depending on the WD mass and the accretion rate 



Nomoto fc Kondo I l|l99ll ) pointed out that, if the accretion 



rate M wd > 10 -8 M yr _1 , an ONe WD with an initial mass 
of 1.2 Mq is very likely to collapse to be an NS rather than 
experience type la supernova explosion. By 2.5-dimensional 
radiation- hydr o dyna mics simulations of the AIC of WDs, 
iDessart et all (|2006l ) proposed that the spin period of the 
natal NSs are 2.2 - 6.3 ms. Therefore, in our calculations 
all NSs formed by AIC are assumed to be rapidly rotating 
objects with millisecond period, and low magnetic field. In 
the input parameters, the initial magnetic field strengths of 
the MSP are chosen from a lognormal distribution of mean 9 
and standard deviation 0.4, and its initial spin periods Pi are 
distributed uniformly between 1 and 10 ms. If the donor star 
fills its Roche lobe after the AIC, the NS should appear as an 
X-ray pulsar (|Stairs 1120041 ; iFreire | [2009). and radio emission 
is suppressed. However, even if for a detached binary with a 
He star or a massive companion, accretion may occur due to 
capture of the stellar winds, which may strongly influence 
the radio emission of MSPs. Base d on the simplified versio n 
of the theoretical model given bv lDavies Sz Pringle I |l98ll ). 
we consider the MSP's spin evolution as follows. 

(i) Radio pulsar phase: After the MSP birth via AIC, it 
appears as a radio pulsar due to the rapid rotation, which 
causes its radiation to be strong enough to expel the winds 
coming from the He star beyond the radius of the light cylin - 
der n c or the Bondi accretion radius rarr llDai et al.l 12006') . 
As a result of magnetic dipole radiation, the spin angular 
momentum loss rate of the MSP is 



Jul — 



2/i 2 fi 3 
3c 3 



(4) 



where fl, and /i = 10 30 ^i3oGcm 3 are the angular velocity, 
and the magnetic dipole moment of the MSP, respectively. 

Considering the interaction between the magnetic field 
of the MSP and t he wi nd material from the donor star, 
iDavies fc Pringlel |l981) proposed that, if the spin period 
P of the NS is greater than either 



1.2M- 1/4 ^V /2 s, 



P ab = o.8Mr 6 1/6 /4o 3 « 8 



5/6 



M 



1/3 



(•») 



(6) 



the radio phase stops. Here M — 10 15 Misgs _1 is the rate of 
mass flow onto the NS, v = 10 8 «gcms _1 is the wind velocity 
relative to the NS, and M is the mass of the MSP. 

Based on the criterion for the radio phase mentioned 
above, a NS formed via AIC is assumed to be a MSP until 
its spin period P > 10ms, or it crosses the so-called "deat h 
line", i.e. B 12 /P 2 < 0.17Gs" 2 l|Bhattacharva et al. |[l993 ). 
where B = 10 12 £>i2G is the surface magnetic field of the NS. 

(ii) Propeller phase: With the spin-down of the NS, the 
light cylinder move outwards, and the magnetosphere radius 



Model 


<*CE 




CAIC 


A 


3 


265 


150 


B 


1 


265 


150 


C 


3 


265 


100 


D 


1 


265 


100 


E 


3 


265 


50 


F 


1 


265 


50 



is less than r\ c , and is greater than the corotation radius 
r c = 1.5 x 10 8 (M/M Q ) 1/3 P 2/3 cm. Subjecting to the cen- 
trifugal barrier, the winds material is assumed to be ejected 
at r m , and exerting a prope ller spin-down torque on the NS 
l|lllarinov fc Sunvaevl ll9T5h . The spin angular momentum 
loss rate via the propeller effect is given by 



j p = 2Mri(l K (r m )[l - Q/n K (r m )], 



(8) 



r m = 1.6 x 10 8 P^ 7 M 18 



2/7 



(7) 



where f2K(V m ) is the Keplerian angular velocity at r m - 

(iii) Accretion phase: we neglect the influence of the 
winds accretion on the spin of the MSP owing to small spin 
angular momentum of wind material. 



3 RESULTS 

Based on the theoretical model described in section 2, we cal- 
culated the evolution of n 3 binaries to an age of 12 Gyr using 
the BSE code, as well as additions for MSP evolution which 
will affect the lifetime, St, used to calculate BMSP numbers 
(see Eq. 3). To explore the effect of input parameters, we 
have calculated the evolutionary results for six models (see 
Table 1), which is determined by the CE parameter oce, the 
kick dispersion tree of core collapse, and the kick dispersion 
caic via AIC channel. We collect the evolutionary stages of 
MSPs as follows: (1) Binary MSPs with a main sequence, 
red giant, or subgiant companion (BMSP-MG); (2) Binary 
MSPs with a He star companion (BMSP-He); (3) Isolated 
MSPs. In Table 2, we summarize the predicted numbers and 
birthrates of various types of MSPs via AIC channel in the 
Galaxy. 

To evaluate the uncertainties in the birthrates and num- 
bers for various types of BMSPs, in table 2 we also present 
the simulated results for same model under different random 
number seeds and grid points n x . We find that, the pre- 
dicted numbers have considerable uncertainties, especially 
the numbers of BMSP-He in models D and F for different 
random number seeds, as well as the numbers of BMSP-MG 
in models A, C, and E for different grid points. However, 
the birthrates besides BMSP-MG show only small varia- 
tions (with uncertainties < 20%). The uncertainties in the 
numbers seem to originate from various random number se- 
quences under different seeds or grid points. The reasons are 
as follows. Firstly, the random number sequence can influ- 
ence the kick received by the natal MSP, which can affect 
the evolutionary fate of binaries. Secondly, the spin evolution 
of MSPs is related to the distribution of the initial period 
and magnetic field. Both parameters strongly depend on the 
random number sequence, which can result in the different 
evolutionary histories of MSPs (see Eq. 3). 
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In model A, we take a strong kick of ctaic = 150 kms 
for the AIC channel, and oce = 3. Our results show that 
there exist ~ 360 BMSPs, which formed via AIC channel. 
In our standard model, BMSPs with a He star companion 
have a number that is ~ 10 — 100 times higher than that 
with a main sequence, red giant, or subgiant companion, and 
the birthrate also shows the same tendency 0. In addition, 
the AIC channel can result in the formation of ~ 5.0 x 
10 s isolated MSP s, which is com patible with the estimation 
by observations (|Lorimenli"995h . As a result of large kick, 
~ 100 - 200 BMSPs with an eccentric orbit (e > 0.1) is 
produced. 

The results of model B show that, for the same kick, a 
lower oce results in fewer BMSPs. Other models also have 
the same tendency. This difference originates from the influ- 
ence of acs on the evolution of close binaries. A high cxce 
can prevent coalescence during the CE stage, significa ntly 
increasing the formation rate of BMSPs (|Liu fc Li l200d ). In 
model E and F, we take a low kick for AIC channel, and 
find that more BMSPs are produced, whereas almost no ec- 
centric BMSPs are formed. In addition, most models predict 
that the number of BMSPs with a RG, or SG companion (no 
BMSPs with MS donor star) is less than 10 (besides model 
A, C, and E when n x = 120), and their eccentricities close 
to 0. 

In Figure 1, we plot the distribution of BMSP-He bina- 
ries in the e — P OI b diagram for models A, B, and C under 
seed 1 and n x = 100. It is seen that in model B most BMSP- 
He binaries are in tight orbits (P or b < Id) with eccentricities 
e > 0.4. For models A and C, the binary population seems 
to be distributed into two regions. Similar to model B, the 
first group have a short orbital period of < Id and a large 
eccentricities of e > 0.4. While contrary to the first group, 
the second are consists of BMSPs with a long orbital period 
of > 20d and an eccentricities near 0. 

The distribution of the primordial binary systems form- 
ing BMSPs in the Mi - M 2 (left) and Mi - a (right) di- 
agrams are shown in Figure 2 for models A, and E. One 
can see that the parameter space of the primordial bina- 
ries forming BMSP-He is wider than that of BMSP-MG, 
and the latter has the rigorous initial parameters. In addi- 
tion, the initial parameter space of BMSP-He in model A is 
larger than that in model E, while the tendency is reversed 
for BMSP-MG. This difference results from the various ef- 
fects of the kicks on BMSP-He and BMSP-MG. The former 
usually has a narrow orbit, and high kick can prevent the 
merger of two components. However, for BMSP-MG with a 
wide orbit high kick can lead to more binaries to be broken 
up during the AIC. The difference in the initial parameter 
space naturally account for the difference of birthrates. BM- 
SPs with a long orbital period (> 20d) have evolved from 
primordial binaries with Mi ~ 8 - 11 M©, M 2 ~ 2 - 4 M Q , 
and a ~ 100 - 1000 R@. 



3 iTutukov fc Yungelson I il996h suggested that binaries with non- 
degenerate He star may be the origin of the most vast AIC, and 
outnumber binaries with H main sequence star. 



4 DISCUSSION AND SUMMARY 

The main goal of this work is to explore if the AIC channel 
can form a population of eccentric BMSPs. For a massive 
ONe WD that accretes material from its donor star, when 
its mass reaches Chandrasekhar limit, the process of electron 
capture may induce gravitational collapse of the WD rather 
than a type la supernova explosion. As a result of AIC, an 
MSP with a weak surface magnetic field, and low spin-down 
rate may be formed. Same as the formation of those NSs 
originated from the collapse of massive star, the nascent 
NSs via AIC may also receive a kick due to the asymmetric 
ejection of material or neutrinos from the newborn NSs. Our 
simulated results can be summarized as follows. 

1. For a high kick velocity (150kms _1 ) during the AIC, 
there exist ~ 100 - 200 (q C e = 3) or ~ 10 - 30 (q C e = 1) 
eccentric BMSPs in the Galaxy, accompanied by a He star. 
However, almost all BMSPs are in circular orbits when the 
natal kick velocity < 100km s -1 . 

2. The ONe WD + He star channel play an important 
role in forming BMSPs. Except for model D, our predicted 
number and birthrate of BMSP-He are 1-2 orders of mag- 
nitude higher than those of BMSP- MG. In fact, this evolu- 
tionar y channel h as been notice d by Tutukov fc Yungelson I 
ljl996t l. Recently, IWang et alD |2009al lbl) investigated the 
evolutionary channel of CO WD + He star to type la su- 
pernovae, and found that this channel can interpret type la 
supernovae with short delay times. 

3. It is difficult for the AIC channel to produce BMSPs 
with a MS, RG, or SG companion, the number of which is 
less than 10, and their eccentricities are near 0. If the com- 
panion of PSR J1903+0327 is really a MS star, the prob- 
ability that it originated from the current binary can be 
ruled out. Both the recycled model and the AIC channel 
cannot produce this peculiar BMSP with a high eccentric- 
ity and a MS donor star, unle ss the natal NS had possessed 
a debris disk (|Liu fc Lil f2009'). We expect further detailed 
multi-waveband observations for PSR J1903+0327 to con- 
strain the properties of its optical counterpart, and hence 
its formation channel. 

4. All the models can produce ~ 10 5 isolated MSPs in 
the Galaxy, without invoking the evaporatio n of donor stars 
by the high energy radiation of the MSPs (Kluznia k et al.l 
ll98StlRuderman et al. Ill989al lbl). In addition, AIC of WDs 
may form sub-m illisecond pulsars in a quark star regime 
jDuet al. Il2009h . which can hardly be formed in the sce- 
nario of normal NSs. In particular, the AIC channel seems 
to be related to the discrepancy between the birthrates of 
MSPs and LMXBs, the solution of which may need addi- 
tional mechan isms forming MSPs (|Bailvn fc Grindlav Ill990l ; 
lLorimerlll995h . Certainly, it is interesting to investigate fur- 
ther in this alternative model. 

5. We note that the predicted numbers of various types 
of BMSPs have considerable uncertainties for different ran- 
dom number sequences. However, the birthrates and the 
numbers for BMSP-He in a given model show little or no 
variation. Since BMSP-He make a significant contribution 
to the total numbers of BMSPs, our results are reliable in 
this respect. 

Perhaps the biggest uncertainty in this work is the kick 
distribution during the AIC. It strongly depend on the kick 
velocity if the AIC channel can form the eccentric BMSPs. It 
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Figure 1. Orbital period and eccentricity probability distribution of BMSPs formed via the AIC channel for models A, B, and C, from 
left to right, respectively. 

Table 2. Predicted numbers and birthrates of various types of MSPs via the AIC channel for differential models in the 
Galaxy. 
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A 


scedl,100 


363 
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1 


4.3 X 10" 6 


5.0 X 10 5 
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seed2,100 
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4.0 X 10 5 
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2.1 x 10 5 
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seed2,100 
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2.0 x 10" 6 


2.6 x 10 5 
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seedl,120 
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5.1 x 10~ 5 
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2.2 x 10 5 
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seedl,100 


301 


2.1 x 10~ 4 
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4.6 x 10~ 6 


4.0 x 10 5 
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seed2,100 
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4.2 x 10" 6 


3.6 x 10 5 


1.8 x 10~ 4 
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sccdl,120 
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3.9 x 10 s 


1.8 x 10" 4 
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seedl,100 
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seed2,100 
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is usually assumed that the natal NS formed by AIC may re- 
ceive a small kick JPodsiadlowski et alJl2004l: IScheck et al. I 



120041 ; iDessart et al. I l200fj ; iKitaura et al.ll2006l ~ Therefore, 
the possibility forming the eccentric BMSPs via the AIC 
channel can be ruled out. In addition, even if the natal NS 



receives a strong kickQ, AIC channel can also not produce 



4 The natal NS via the AIC would los e ~ O.2M0 mass in the form 
of neutrinos jvan Para diis et al. 1997). If the ejection of neutrinos 
is highly asymmetric, a large kick may be produced. Alternatively, 
during AIC a large amount of mass ejection may occur as a result 
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Figure 2. Distribution of primordial binary systems forming BMSPs in the Mi — M2 and Mi — a diagrams for models A (top panels), 
and E (bottom panels). The open squares, crosses, and solid circles correspond to the progenitors of BMSP-MG, BMSP-Hc, and BMSPs 
with a long orbital period, respectively. 
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the eccentric BMSPs with an orbital period of > 10 days, 
and their companions are He stars. Based on the results by 
binary population synthesis, we demonstrate that the ob- 
served orbital parameters of PSR J1903+0327 cannot be 
produced by the AIC channel. 
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